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^SJ , Abstract 

^ ■ We study the classical and quantum models of a Friedmann-Robertson- Walker (FRW) cosmol- 

^■-N I ogy, coupled to a perfect fluid, in the context of the f{R) gravity. Using the Schutz' representation 

■ for the perfect fluid, we show that, under a particular gauge choice, it may lead to the identification 

of a time-parameter for the corresponding dynamical system. Moreover, this formalism gives rise 
to a Schrodinger-Wheeler-DeWitt (SWD) equation for the quantum-mechanical description of the 
I model under consideration, the eigenfunctions of which can be used to construct the wavefunction 

of the Universe. In the case of f{R) — (pure quadratic model), for some particular choices of 
' the perfect fluid source, exact solutions to the SWD equation can be obtained and the correspond- 

ing results are compared to the usual f{R) — R model. 

CN ■ PACS numbers: 98.80.Qc, 04.50.-t-h, 04.60.Ds 
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Over the past few years, extended theories of gravity, constructed by adding higher-order curvature 
• terms to the usual Einstein-Hilbert action, have opened a new window to study the accelerated 

expansion of the universe. It has been shown that such correction terms could give rise to accelerating 
solutions of the field equations without having to invoke concepts such as dark energy [1] . In a more 
general setting, one can use a generic function f{R), instead of the usual Ricci scalar R, as the action 
^ ■ of gravitational field. In more recent times, such f{R) gravity theories have been extensively studied 

in the literature, see [2] for a review. 

In fact, quadratic and higher-order gravity theories have been used long before the accelerated 
expansion of the Universe is observed, immediately after Lovelock's natural generalization of the 
Einstein-Hilbert action for the gravitational field [3]. Quadratic Lagrangians, in particular, have 
been used to yield renormalizable theories of gravity coupled to matter [4]. They can also help us 
to improve the semiclassical approximation, where quantized matter fields interact with the classical 
gravitational field [5]. In this context, renormalization of the energy-momentum tensor for a quantum 
field in a four-dimensional, curved space-time indicates that the presence of quadratic terms in the 
gravitational action is a priori expected [6]. On the other hand, Kleidis et al [7], demonstrated that, 
in quadratic gravity theories, an additional coupling arises between i?^ and any massive quantum 
scalar field, introducing a geometric source term in the wave equation for the quantum field. 

In order to find the dynamical equations of motion one can vary the action with respect to the 
metric (metric formalism), or view the metric and connections as independent dynamical variables 
and vary the action with respect to both independently (Palatini formalism) [8]. One should note 
that, although in the usual Einstein-Hilbert action these two approaches give the same field equations, 
in f{R) gravity the Palatini formalism leads to different dynamical equations due to nonlinear terms 
in the action. However, the f{R) form of the action is shown to be equivalent to a scalar-tensor type 
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theory and there is an equivalence between the metric and the Palatini f{R) gravity with Brans- 
Dicke theories. Indeed, one can show that these versions of f{R) gravity are dynamically equivalent 
to the Brans-Dicke theory of gravity with values of the Brans-Dicke parameter u> = and uj = —3/2, 
respectively [2], [9]. There is also a third version of f{R) gravity in which the Lagrangian of the 
matter depends on the connections of the metric (metric-affine formalism) [10]. 

It becomes evident that, the corrections of the /(i2)-gravity to the results of the standard general 
relativity have been widely investigated, but, most of these works are restricted in the classical regime 
[11]. The cases dealing with quantum f{R) models have seldom been studied in the literature [12], 
and it would be of interest to employ such models in this study. 

In this paper we shall investigate the problem of quantization of a FRW cosmological model in 
the context of f{R) cosmology. For the matter field coupled to gravity we consider a perfect fluid in 
Schutz' formalism [13]. The advantage of using this formalism in our quantum cosmological model is 
that, in a natural way, it can offer a time parameter in terms of dynamical variables of the perfect fluid 
[14]. Indeed, as we shall show, after a canonical transformation the conjugate momentum associated to 
one of the variables of the fluid appears linearly in the Hamiltonian of the model. Therefore, canonical 
quantization results in a Schrodinger-Wheeler-DeWitt (SWD) equation, in which this matter variable 
plays the role of time. However, it should be emphasized that we do not deal with the problem of time 
in the quantum cosmology of the model at hand in a fundamental way. As are well known, quantum 
cosmology suffers from a number of problems which one of the most important of them is the problem 
of time, i.e., the wave function in the Wheeler-DeWitt equation is independent of time and thus the 
Universe has a static picture in this scenario. This problem was first addressed in [15] by DcWitt 
himself. He argued that the problem of time should not be considered as a hindrance in the sense that 
the theory itself must include a suitable well-defined time in terms of its geometry or matter fields. 
Therefore, since by the Whcclcr-DeWitt approach one can find different time parameters directly 
from the the super-metric, the sole purpose of the use of a perfect fluid instead of a field with its own 
degrees of freedom seems to be the emergence of a Schrodinger-type equation with a time parameter 
t, which is not the case for the general Wheeler-DeWitt equation. The Schutz model is thus of limited 
value in view of the problem of time in quantum cosmology. 

Here, we first consider a generic form of f{R) function in the FRW background with a generic 
equation of state for the perfect fluid to construct the phase space of the model. As we find out, even 
for the simplest forms of the /(i?)-function, the corresponding SWD equation is too complicated to 
be solved in terms of analytic functions. Therefore, we restrict ourselves to a pure fourth order model, 
in which f{R) = R^ ■ Accordingly, we show that, for some special forms of the equation of state of 
the perfect fluid, the SWD equation can be solved exactly in terms of known special functions. We 
then construct the resulting wavefunction of the universe in this quadratic framework and compare its 
similarities and differences with usual f{R) = R model. It is to be noted that our presentation does 
not claim to deal with f{R) quantum cosmology in a fundamental way since we study the problem in 
a simple toy model. However, this study may reflect realistic scenarios in other similar investigations, 
dealing with such problems in a more fundamental way. 

2 The phase space of the model 

In this section we consider a FRW cosmology within the framework of f{R) gravity. In the context 
of the ADM formalism, the action for gravity coupled to a perfect fluid in Schutz' representation is 
written as (in what follows we work in units where c = h= IGttG = 1) 

S= f (fx^f{R) + 2 [ d^xVhhabK"^ + I (fx^p, (1) 

JM JdM JM 

where R is the scalar curvature and f{R) is an arbitrary function of R. Also, is the extrinsic 
curvature and hab is the induced metric over the three dimensional spatial hypersurface, which is 
the boundary dM of the four dimensional manifold M. The last term of (1) denotes the matter 
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contribution to the total action where p is the pressure of perfect fluid which is hnked to its energy 
density by the equation of state 

P = ap. (2) 

In Schutz' formahsm [13], the fluid's four- velocity is expressed in terms of five potentials e, C,, (3, 9 
and S as 

= -(e,. + CA- + 0S,), (3) 

where ijl is the specific enthalpy, the variable S is the specific entropy while the potentials C and (3 are 
related to torsion and are absent in the FRW models. The variables e and 9 have no clear physical 
interpretation in this formalism. The four-velocity satisfies the condition 

U^U^, = -1. (4) 

We assume that the geometry of spacetime is described by the FRW metric 



ds' = -N\t)dt' + a\t) 



1 - kr' 



+ r^(dr + sin^ 



(5) 



where Nil:) is the lapse function, a{t) the scale factor and A;=l, and —1 correspond to the closed, 
flat and open universe respectively. To proceed further, we need an eS'ective Lagrangian for the model 
whose variation with respect to its dynamical variables yields the appropriate equations of motion. 
Therefore, by considering the above action as representing a dynamical system in which the scale 
factor a, scalar curvature R and fluid's potentials play the role of independent dynamical variables, 
we can rewrite the gravitational part of action (1) as 



Sgrav = j dtCgrav{a,d,R,R) = j dt I Na^f{R) - A 



6 /a k Na^ 



(6) 



where we have introduced the definition of R in terms of a and its derivatives as a constraint. This 
procedure allows us to remove the second order derivatives from action (6) . The Lagrange multiplier 

A can be obtained by variation with respect to R, that is, A = Na"^ f (R), in which a prime denotes the 
derivative with respect to R . Thus, we obtain the following point-like Lagrangian for the gravitational 
part of the model 



J^grav = " ^a'a^/" + 6fc7Va/' + Na^f - Rf). (7) 



To simplify this Lagrangian, we define the variable (p as f'{R) = (p, in terms of which the Lagrangian 
(7) reads 

Cg^^^ = -—aa^cf) - j^a^a^ + QkNa(f) - Na^V{(t)), (8) 

where V{(j)) = Rf - f = Rcj) - f . 

Also, the matter part of the action (1) becomes Smatt = I d^xdtNa^p, so the Lagrangian density of 
the fluid is Cmatt = Na^p. Following the thermodynamic description of [13], the basic thermodynamic 
relations take the form 

p = Po(i + n), fi = i + u + ^, (9) 

Po 

where po and 11 are the rest-mass density and the specific internal energy of the fluid respectively. 
These quantities together with the temperature of the system r, obey the first law of the thermody- 
namics rdS = dll + pd(l/po), where can be rewritten as 



rdS = dU + pd 



(-) =(l + n)d[lnpo"(l + n)l, (10) 
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in which we have used the equation of state (2). Therefore, we obtain the fohowing expressions for 
the temperature and the entropy of the fluid 



r = l + n, 5 = lnpo"(l + n) =ln 



a + 1 



Now, we can express po ^.nd 11 as functions of /x and S as 



i + n = 



Po 



a + 1 



l/a 



-S/a 



a + 1' 

so that with the help of (9), one can put the equation of state in the form 

a 



P 



(11) 



(12) 



(13) 



On the other hand, normalization of the fluid's four- velocity (3), according to the relation (4) implies 
// = (e + 6S)/N. Therefore, using the above constraints and thermodynamical considerations for the 
fluid we find 



matt 



a 



e + es 



l+l/a 



5/a 



(14) 



(a + l)i+Vo 

Let us now construct the Hamiltonian for our model. The momenta conjugate to each of the above 

3 _ ac 

1 ~ dq- 



variables can be obtained from the definition Pn = In terms of the conjugate momenta the 



Hamiltonian is given by 



H = H. 



grav 



+ H, 



matt 



aPa + cl)Ps + ePe + SPs - jC, 



(15) 



where £ = Cgrav + ^matt- Noting that 



P.= 



de 



= N-y 



l/a 



(a + 1)V« 



-s/a 



and 



dS 



(a + 1)V« 



expression (15) leads 

H = Nn = N 



PaP<p 



+ 



6a^ 6a^ 



(16) 



Now, consider the following canonical transformation which is a generalization of the ones used in 
[16] 



T = -Pse-^Pe 



-(a+1) 



(17) 



e = e - (a + 1)-^, 



P. 



Under this transformation Hamiltonian (16) takes the form 



H = NH = N 



6a2 



+ -^P^-^kacl) + a'V{4,) + 



,3a 



(18) 



We see that the momentum Pt is the only remaining canonical variable associated with matter and 
appears linearly in the Hamiltonian. The setup for constructing the phase space and writing the 
Lagrangian and Hamiltonian of the model is now complete. 
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The classical dynamics is governed by the Hamiltonian equations, that is 



f a = {a,H} = -f^, 

Pa = {Pa,H} = N f-|a-3p„P^ + i0a-4p| + 6k(P - 3a'^V{(P) + Saa-^^'-^PT 



<P = {cl),H}=N 
P^ = {P^,H} = N 



1 Pa I 1 </> ; 

'6^ + 3^^ 



(19) 



iJ + 6ka-a^V'{<P) 



{T,H} 



N 



[ Pt = {Pt,H} = 0. 

We also have the constraint equation 7-1 = 0. Up to this point the cosmological model, in view of the 
concerning issue of time, has been of course under-determined. Before trying to solve these equations 
we must decide on a choice of time in the theory. The under-determinacy problem at the classical level 
may be removed by using the gauge freedom via fixing the gauge. A glance at the above equations 
shows that choosing the gauge N = a^", we have 



(20) 



which means that variable T may play the role of time in the model. Therefore, the classical equations 
of motion can be rewritten in the gauge A'^ = a^" as follows 

a = -ia3"-2p., 



P„ = _ la3a-3p^p^ ^ 1^^3a-4p2 ^ QJ^^^a^ _ Sa^'^+^V {(f)) + Saa'^Po, 
<j) = - ia3«-2p„ + la3"-3</,P^, 
P^ 



(21) 



where we take Pt = Pq = const, from the last equation of (19). Since integrability of this system 
directly depends on the choice of a form for f{R) (which determines potential V{4>)), it is appropriate 
to concentrate on this point first. However, before choosing such a function, let us deal with the 
quantum cosmology of the model described above. 



3 Quantization of the model 

We now focus attention on the study of the quantum cosmology of the model described above. We 
start by writing the Wheeler-DeWitt equation from Hamiltonian (18)^. Since the lapse function N 
appears as a Lagrange multiplier in this Hamiltonian, we have the Hamiltonian constraint H = 0. 

^The canonical transformation (17) is applied to the classical Hamiltonian (16), resulting in Hamiltonian (18) which 
we are going to quantize. To make this acceptable, one should show that in the quantum theory the two Hamiltonians 
are connected by some iniitary transformation, i.e. the transformation (17) is also a quantum canonical transformation. 
A quantum canonical transformation is defined as a change of the phase space variables {q,p) — » {q' ,p') which preserves 
the Dirac bracket [17] 

[q,p] = i = [q'{q,p),p'{q,p)] . 
Such a transformation is implemented by a function C{q,p) such that 

q'{q,p) = CqC~^, p {q,p) = CpC~^ . 
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Thus, application of the Dirac quantization procedure demands that the quantum states of the 
Universe should be annihilated by the operator version of 7Y, that is 



H*(a,0,r) 



*(a,(/.,T) = 0, 



(22) 



where ^'(a, ^, T) is the wavefunction of the universe. A remark about the reduced Hamiltonian in 
the above procedure is the factor-ordering problem when one embarks on constructing a quantum 
mechanical operator equation. In dealing with STich Hamiltonians at the quantum level extra care 
should be taken when one tries to replace the dynamical variables with their quantum operator coun- 
terparts, that is , in replacing a variable q and its momentum pg with their corresponding operators, 
the ordering considerations should be taken into account. Therefore, to guarantee Hermiticity, the 
operator form corresponding to equation (22) may be written as 



-6ka(t) + a^y(0) + ^ 



+ <p''"P4>r 

^ia,cf>,T) = 0, 



(23) 



where the parameters r, s, u, v and w satisfy r + s = —2, u + v + w = 1 and denote the ambiguity 
in the ordering of factors a and Pa in the first and 4> and in the second term of (22). With the 
replacement Pa ^ —i-^ and similarly for and Pt the above equation reads 



1 

6"" 



d 



dadcf) 



1 



1 



d<j) 6"' "'a^^ 

+a'V(0) - a-'^'-H— 



6 

d_ 
df 



Gkaq 



*(a,</),r) = 0, 



(24) 



where A = uw. This equation takes the form of a Schrodinger equation id^/dT = H^, in which the 
new Hamiltonian operator is Hermitian for any choice of the ordering parameters with the standard 
inner product 

< # I >= f a-^''^*^dad(t>. (25) 
We separate the variables in the SWD equation (24) as 

*(a,0,r) = e-^^^V(a,</'), 

leading to 



(26) 



a(f) 



^2 



dadcj) 



A - 36A;a^r + 6a>y((^) - 6Ea''-'"^(l) } ip{a, 0) = 0, 



(27) 



where £' is a separation constant. It is seen that this equation has a mixed derivative with respect 
to the variables a and ^ and also, these variables appear in a mixed form in the last term of the 
equation. Under these conditions equation (27) cannot be solved by the method of separation of 
variables. Therefore, it is useful to introduce the following change of variables 



(28) 



This canonical transformation C transforms the Hamiltonian as H'{q,p) = CH{q,p)C ^. For our case the canonical 
relations [S, Ps] = \e, Pc] = i yield 



[T,Pt1 = [-Pse-'Pr'~"+'\P:+'e'] = - [Pse'^e^] = [e^Ps] e 



= le e 



and 



[e, P.] = [e- (a + l)PsP-\P,] = [e, P,] = i, 

which means that the transformation (17) preserves the Dirac brackets and thus is a quantum canonical transformation. 
Therefore, use of the transformed Hamiltonian (18) for quantization of the model is quite reasonable. 
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In terms of these variables equation (27) takes the form 



1 



92 



1 d 



92 



d 



4 dx^ 4 dx 



dy 



-x^— - -X— - j/2 - 2y— + A- 36kx^ + 6x^y" - 6£;x-^--^"y 2 (3"-i) I ^(x, y) = 0. (29) 



Unfortunately, this equation cannot be solved analytically for an arbitrary potential function V{y) 
and parameter a, which represent the form of f{R) and the nature of the perfect fluid respectively. 
In the next section we shall present a class of exact solutions for this equation in pure cosmology 
for two special cases of perfect fluid. 

4 Quantum solutions for f{R) = cosmology 



For the pure fourth order {f{R) = R?), model we have 



f'{R) = 2R and F((/)) = Rf - f = R^. 



Hence, from (28) we obtain F(^) = = iy2_ Therefore, equation (29) reduces to 



1 2^ 

4^ 



1 d_ 



2d' 

y ^ 



5y2 



d 

''e-y^ 



V'(x,y)=0. (30) 



It is seen that in the cases of a = 1, 1/3, the variables x and y can be separated from each other 
and equation (30) may be solved analytically for these two cases. Thus, in what follows we restrict 
ourselves to these two special cases. 



4.1 Stiff matter: a = 1 

Stiff matter is a fluid with pressure equal to energy density and speed of sound equal to speed of light. 
In this case we separate the solutions of equation (30) into the form i/j{x,y) = X{x)Y{y) leading to 



„2 



- + (1 - z/2 - UAkx'^ + 6x^^)1 X{x) = 0, 



(31) 



where we take ^^-j^ as a separation constant and also ^4 = without losing general character of 
the wavefunction. One should note that there exists an infinite number of possibilities of ordering. 
As Hawking and Page have shown [18], the factor ordering parameter will not affect semiclassical 
calculations in quantum cosmology and so for convenience one usually choose a special value for it 
in the special models. On the other hand, in general, the behavior of the wavefunction depends on 
the chosen factor ordering [19]. But in the model at hand, as is clear from equation (30), the factor 
A will appear in the solutions of (31) only together with the separation constant v (for example 
V — V — 4^ in the first equation of (31)). Therefore, since we finally take a superposition over 
all values of this separation constant, we expect that the essential features of the wavefunction will 
remain the same for different choices of the factor ordering parameter A. Also, with an eye to equation 
(23), it is to be noted that the choice ^4 = yields a zero value for u or w; or both of them. In all 
of these cases the Hermiticity of the second term in equation (23) will remain unchanged. Prom this 
point of view the particular choice A = does not reflect major effects to equation (23). 

Por the flat PRW metric {k = 0), the above equations have the following solutions in terms of 
Bessel functions 



X{x) = X [ci J| (y|x3) + csYi (yfx^) 



. ^(y) = ^ [diJ>^ iV^) + d2Y. (x/24E^] , 



(32) 
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where Cj and {i = 1, 2) are integration constants. Thus, the eigenfunctions of the SWD equation 
for stiff matter, can be written as 

*,B(x,y,T) = (V?""') (^^^) ' ^^^^ 

where we choose C2 = d2 = for having well-defined functions in all ranges of variables x and y. We 
may now write the general solutions to the SWD equations as a superposition of the eigenfunctions, 
that is 

POO POO 

^{x,y,T)= / A{E)C{u)^,E{x,y,T)dEdu, (34) 

Je=o Ju=0 

in which A{E) and C(i^) are suitable weight functions to construct the wave packets. By using the 
equality [20] 

70 ^ (2a)-+i' 

we can evaluate the integral over E in (34) and simple analytical expression for this integral is found if 
we choose the function A{E) to be a quasi-Gaussian weig ht factor A{E) = 12{2AEY/^e-^'^"'^ , which 
results in 

where 7 is an arbitrary positive constant. Substitution of the above relation into equation (34) leads 
to the following expression for the wavefunction 

^ix,y,T) = ^exp (-—^] H , ^^'l.u^i J^. { xflA (37) 



00 „ . 1,2 



g-ar- j.^+^j^{i)j.)dr = \ e-fe, (35) 



(27 + S) 



where C{v) can now be chosen as a shifted Gaussian weig ht function e-"(''-^)'. It is seen that this 
expression is too complicated for extracting an analytical closed form for the wavefunction and the 
choice of a function C{i') that leads to an analytical solution for the wavefunction is not an easy 
task. In this respect, our choices for A{E) and C{v) as quasi-Gaussian and shifted Gaussian weight 
functions respectively, appear to have physical grounds. Such Gaussian weight functions are widely 
used in quantum mechanics as a way to think about the localized states. This is because that these 
types of weight factors are centered about a special value of their argument and they fall off rapidly 
away from that center. Due to this behavior the corresponding wave packet resulting from (37) 
after integration, has also a Gaussian- like behavior, i.e., is localized about some special values of its 
arguments. Therefore, it is seen that there is a reciprocal relation between the width of the Gaussian 
weight function that determines the shape of the wave packet, and the width of the wave packet. In 
this sense the motion of the peaks of the wave packets with a group velocity and also its spreading 
behavior (see below), can be best seen in terms of the Gaussian packets. 

In figure 1 we have plotted the square of the wavefunction for typical numerical values of the 
parameters. As this figure shows, at T = 0, the wavefunction has a dominant peak in the vicinity 
of some nonzero values of x and y followed by smaller peaks, which as x grows, their amplitudes are 
suppressed. This means that the wavefunction predicts the emergence of the universe from a state 
corresponding to its dominant peak. However, the emergence of several peaks in the wave packet may 
be interpreted as a representation of different quantum states that may communicate with each other 
through tunneling. This means that there are different possible universes (states) from which our 
present universe could have evolved and tunneled in the past, from one universe (state) to another. 
As time progresses, the wave packet begins to propagate in the y-direction, its width becoming wider 
and its peaks moving with a group velocity towards the greater values of y. Bearing in mind that 
y = 4> = f'{R)-, tliis wavefunction predicts that the universe will assume states with larger R in its 
late time evolution. 
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2.5 2.5 



Figure 1: The figures show \'if{x,y,T)\^, the square of the wavefunction in four different time parameter T — 
0, 10, 15, 20. The figures are plotted for the numerical values a — 2, b — 5 and 7 — 0.5, and we have taken the 
integral in (37) from to 10. 
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Figure 2: Approximate behavior of a{t), the scale factor, x{t) = atf)^^^ and R{t), the Ricci scalar, in the classical 
cosmology framework. The figures are plotted for numerical values a — 1, k — and Po = 1. We take the initial 
conditions a{t = 0) = (j>{t = 0) = 1 and Pa{t = 0) — P^{t — 0) = 3. After examining other sets of initial conditions, we 
verify that this behavior repeats itself. 



Before going any further, let us take a look at the classical cosmology corresponding to the above 
quantum solutions. For the quadratic model, we insert V{(j)) = jcj)"^ in the system (21), take a = 1 for 
the matter character and also k = representing the flat FRW space time. We see that the classical 
cosmology forms a system of nonlinear coupled differential equations which unfortunately cannot 
be solved analytically. In figure 2, employing numerical methods, we have shown the approximate 
behavior of a(t), x{t) and (p{t) ^ R{t) for typical values of the parameters and initial conditions 
respectively. As is clear from the figures, while the Ricci scalar has a monotonically growing behavior 
with time, the scale factor first decreases and then increases forever. On the other hand, take a look 
at the figure 1 shows that the quantum states with larger (smaller) peaks, i.e., the states of high 
(low) probability, are located at the smaller (larger) values for the scale factor. Therefore, when the 
scale factor decreases (increases) with time in the classical point of view, in the quantum domain 
tunneling occurs from a state with a smaller (larger) peak to another with a larger (smaller) one. 
Also, in view of the behavior of the Ricci scalar, the classical and quantum solutions are in complete 
agreement with each other and both predict a monotonically increasing evolution for this variable. 
Therefore, it is seen that there is an almost good correlation between the quantum patterns shown in 
figure 1 and classical trajectories. However, in order to show this correlation in a more clear sense, 
one may calculate the time dependence of the expectation values. As we have seen in the second 
section, the variable T can be regarded as an internal clock for the dynamics. The effective (reduced) 
Hamiltonian for this system can be found by solving the constraint TC = for zizPT [21]. Hence 

Heff = -la^'^-^PaP^ + Ja3°-3^p2 _ g^^Sa+l^ ^ a^'^+'^V {(j)) . (38) 

With respect to this Hamiltonian the time evolution of the expectation value of an observable O 
reads as = —i < [0,Heff] >. Therefore, we get 

d<a> _ _1 ^ „3a~2p ^ 
dT ~ 6 ^ ^(p ^1 

-3 < a3"+V(0) > +3aPo < >, 

(39) 



dT ~ 6 

= -1 < a3-3p2 > < „3a+l > _ < a3"+V(0) > . 

In the case of a very localized wave packet as we have constructed in figure 1 (see also figure 3 in the 
next subsection) a function 0(0, P) varies very little over the size of the wave packet and thus the 
approximation < Q(0,P) >~ il(< O >,< P >) holds. Under this condition the trajectories for the 
expectation values obtained from the system (39) are in agreement with the classical ones obtained 
from the system (21) and shown in figures 2, 4. To summarize, in the present case where we have 
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localized wave packets the expectation value of the scale factor is very close to the classical one and 
thus < a > follows the pattern of a{t) shown in figure 2 with a good approximation. 



4.2 Radiation: a — 1/3 



In the standard Big-Bang model, it is believed that the early universe is dominated either by radiation 
or radiation and high energy particles. Since in this era the temperature is above electron-positron 
pair threshold, these particles are relativistic and the whole mixture behaves more like radiation 
rather than matter. When the cosmological fluid is dominated by radiation, as it was presumably the 
case in the early universe, the equation of state can be taken as p = p/3. In this case, the solutions 
of equation (30) may be again separated into the form 'tp{x,y) = X{x)Y{y) leading to 



X 



2 d2 



' dx 



+ (1 - 4i^2 _ 144^2;^ + _ 24£;x2)l X(x) = 0, 



(40) 



2d^ 



y{y) = 0, 



where we take 1 — Av^ as a separation constant. Similar to the previous subsection, we deal again 
with the solutions of the above system in the case when A; = 0. In general, although the second 

equation of the system (40) has exact solutions, its first equation cannot be solved analytically. 
However, as we mentioned above, the radiation is not the dominant fiuid at the present epoch but 
was dominant in early universe, i.e., in the regime where the scale factor is expected to have a small 
value. On the other hand, our classical analysis in this case (see below and figure 4) will show that 
the Ricci scalar follows a vanishing value. Therefore, we may write x"^ = a?{2R) ~ ^ <C 1 and hence 
x^ = a^{2R)^ ^ ^a^(2i?)^ = ^x^. Thus, we use the approximation x^ ~ in this era and rewrite 
the SWD equation as ^ 



X 



,2 d^ 



x^ + (1 - 4z/2 + e^x^ - 2AEx'^) X{x) = 0, 



(41) 



,2 d-^ 



Y{y) = 0. 



Now, the first equation of the system has well known solutions in terms of confiuent hypergeometric 
functions M(a, 5;x) and U{a,b;x). Therefore, we get 



CiU + iVeE, 2v + 1; i^f^xA + C2M f + J - ^V6^, + 1; i^x^yii2) 



X{x) 



Y{y) = D^y^-^l'' + D2y-''-^'\ (43) 

where Q and Di {i = 1,2) are integration constants. Thus, the SWD equation for radiation fluid has 
the following eigenfunctions 



*,^;(x,y,r)=e-^^^x2'^+ie- 



''m{v+]^- iVeE, 2v + 1; i^/^x'^^ y 



-1/-1/2 



(44) 



where we take Ci = Di = to have well-defined eigenfunctions for constructing the wave packets. To 
construct such wave packets since the weight functions can be chosen as a shifted Gaussian function, 
the general solutions to the SWD equation read as 



/•OO /"OO 

^'(x,y,r)= / / , 

Jv=Q JE=0 



-cr(l/-l/o)" 



g-?(E-£o)2^^^(3.^ y, T)dvdE. 



(45) 



In figure 3 we have plotted the square of the above wavefunction for some numerical values of the 



^Notice that this approximation leads to i? ~ which for a constant a and very small value for ^ goes to zero. 
The same result is also achieved from the definition of R in (6) when the scale factor is constant. 
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Figure 3: The figures show j/, r)p, the square of the wavefunction in two different time parameter T — 0, 10. 

For larger T's the same pattern is repeated. The figures are plotted for the numerical values ^ = 0.1, a = <; — 10 and 
vq = Eq = 1, and we have taken the integrals in (45) from to 15 for E and from tO 4 for v. 



-100 ^ 



-20 



y=0~R 



80 




4 



Figure 4: Approximate behavior of a(i), the scale factor and R{t), the Ricci scalar, in the classical cosmology 
framework. The figures are plotted for numerical values a — 1/3, fc = and Po = 1- We take the initial conditions 
a(t = 0) = 1, 0(t = 0) = 0, Pa{t = 0) = —10 and -P^(t = 0) = —1. After examining some other acceptable sets of initial 
conditions, we verify that this behavior is almost repeated. 

parameters. The figures show that the wavefunction has several peaks which are sharply distributed 
around y = 0. Also, it is seen that as time grows, the pattern of the wave packets do not show a 
major change. This means that in the context of our quadratic model the universe in its radiation 
era has almost a constant scale factor and thus its scalar curvature vanishes. 

Now, let us pursue the corresponding classical cosmology given by the system (21) with a = 1/3, 
V{(j)) = \(t? and /c = 0. In fi gure 4 we have prepared the approximate behavior of a{t), and (f){t) ~ R{t) 
for typical values of the parameters and initial conditions respectively. As the figure shows, in complete 
agreement with the quantum model, the universe follows a zero Ricci scalar while its scale factor has 
a constant value in the epoch of the radiation dominated era of cosmic evolution. 

5 Comparison with f{R) = R model 

In this section we briefly review the problem in the usual f{R) = R framework which originally 
appeared in [14]. The Einstein-Hilbert Hamiltonian of the model used in the present work can be 
written as 




(46) 
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Therefore, the equations of motion become 



a = {a,H} = -^, 

Pa = {Pa, H] = -N(^-6k- 3a^) , 

(47) 

f = {T,H} = ^, 
Pt = {Pt,H} = 0. 

Choosing the gauge N = a^" , we have T = t and the following constraint equation H = 

-6-3^ - 6A;a3"+i + Po = 0, (48) 
where Pt = Po = const, as before. For flat space-time (A; = 0), the above equation has the solution 

a{t) = aoe V « , (49) 

for stiff matter (a = 1), and 




a{t) = \^t + ao, (50) 



for radiation (a = 1/3)^. It is clear that for both of these solutions the Ricci scalar R = ^ (^^ + ^ — 
has a behavior R ^ as t ^ +00. In the case when the universe is filled with stiff matter, comparison 
of the solution (49) with those of the quadratic model which appear in figure 2, shows that while 
the scale factors have the same behavior at late times, their early-time evolution does not coincide 
with each other. This means that although the quadratic model can predict the late time expansion 
of the universe, but in spite of the usual f{R) = R model its Ricci scalar grows monotonically with 
time. Also, for radiation the results of two models show some differences. In f{R) = R cosmology 
the scale factor has a linearly expanding behavior according to (50), with a decreasing Ricci scalar 
while the quadratic cosmology shows a static universe with zero Ricci scalar. To give an estimate why 
the result of the quadratic R? -cosmology decline from those of the linear Einstein-Hilbert theory, we 
should note that although this form of f{R) gravity yields a late time expansion cosmology, but in 
view of having the correct weak-field limit at Newtonian and post-Newtonian levels has not a desired 
form. The conditions under which a modified gravity model passes the local and astrophysical tests 
such as Newton law and solar system tests are investigated in [22]. In these works such f{R) theories 
are studied which satisfy the conditions 

lim f{R) = cons., lim /(i?) = 0, (51) 

it— »0O it-+0 

and shown that they pass Newton law, stability of Earth-like gravitational solution, heavy mass for 
additional scalar degree of freedom, etc. Therefore, since our quadratic model does not satisfy the 
above conditions, it is not a viable theory with correct Newtonian and post-Newtonian limits. This 
is not surprising since it is well known that a large class of f{R) theories suffer from this issue [9]. In 
summary, in the framework of quadratic model, although we have a compatible classical and quantum 
cosmology discussed in previous sections, the predictions of this model do not coincide with the usual 
f{R) = R case. 



'in terms of cosmic time f] = J Ndt, these solutions read as a{r]) ~ rj^^^ for stiff matter and a{r}) ~ i]^^^ for radiation. 
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6 Conclusions 



In this paper we have studied a quantum cosmological model in the framework of quadratic f{R) 
gravity coupled to the Schutz' perfect fluid. The use of Schutz' formalism for perfect fluid allowed us 
to obtain a SWD equation in which the only remaining matter degree of freedom plays the role of 
time parameter in the model. For stiff fluid (a = 1) and radiation (a = 1/3), we saw that the SWD 
equation can be separated and its eigenfunctions can be obtained in terms of known special functions. 
In the case of stiff matter as the cosmic fluid, we found the eigenstates in terms of Bessel functions 
and then constructed the corresponding wave packets by appropriate superposition of the eigenstates. 
The wavefunction in this case shows a pattern in which there are different possible quantum states 
(with different probability) from which our present universe could have evolved and tunneled in the 
past from one state to another. The time evolution of this wave packet represents its motion along 
the larger R direction. We have also solved the corresponding classical cosmology in this case by 
numerical methods and showed that the Ricci scalar has a monotonically growing behavior with time 
while the scale factor first decreases and then increases forever. Therefore, we have a model with 
compatible classical and quantum cosmology which although can predict the late time expansion of 
the universe, but contrary to the usual standard model of cosmology the universe will admit states 
with larger R in its late time evolution. On the other hand, in the case of radiation, we obtained 
wavefunctions which preserve their configuration during time evolution. These wavefunctions peak 
sharply around the zero Ricci scalar. The classical cosmology in this case is again in agreement with 
the quantum model and shows an almost static universe with zero scalar curvature. Since in the 
standard model the universe has an expansion behavior as a{r]) ~ 77^/^ (77 is the cosmic time) in its 
radiation dominated era, our quadratic model does not coincide with real cosmological observations 
also in this epoch. 
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